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Abstract 

We explore a class of dark matter models with two dark matter candidates, only one interacts 
with the standard model sector. One of the dark matter is thermalized with the assistance of the 
other stable particle. While both stable particles contribute to the total relic density only one can 
elastically scatter with nuclei, thus effectively reducing the direct detection rate. 
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I. INTRODUCTION 



The most natural explanation for the astrophysical and cosmological indications of a large 
component of dark matter in the Universe is a new weakly interacting massive particle that 
behaves as cold dark matter. Such dark matter candidates are found in different extensions 
of the standard model (SM) However, dark matter (DM) could be composed of more 
than one particle, and several multi component DM models have been suggested recently 
8(1 . Interest for multi component DM arose, in particular, from hints for signals in cosmic 
rays as well as in direct detection experiments corresponding to two completely different 
DM mass scales. On the one hand, the PAMELA experiment 9| has reported an cosmic ray 
excess of positrons which would be consistent with a TeV scale DM candidate annihilating 
mainly into leptons. On the other hand, DAMA 10|, CoGeNT [ll| and CRESST 12| have 
all reported excesses in the direct detection rate that would be compatible with DM around 
10 GeV. Even though there is no conclusive evidence that these signals are due to DM, it is 
interesting to investigate multi component DM models to expand the range of DM models 
at a time when DM searches in indirect, direct and collider experiments are increasing 
their sensitivities. Furthermore, replacing the usual R-parity symmetry that guarantees the 
stability of the lightest R-parity odd particle by an enlarged symmetry group allows not 
only multiple DM candidates but also new freeze-out mechanisms. In particular, the semi- 
annihilation mechanism where two DM particles annihilate into another DM particle and a 
SM particle was proposed in Ref. 13|. 

In this paper, we propose a new type of freeze-out mechanism, assisted freeze-out, within 
the framework of multi component DM models. In this new freeze-out mechanism, one DM 
candidate can be thermalized only through the assistance of the other stable particle. Thus, 
the decoupling of one DM particle from the thermal bath is influenced by the other DM 
particle. Consequently, the relic density of DM is solved by using two coupled Boltzmann 
equations for two stable fields. In the analysis of the right-handed sneutrino DM [l^, a 
similar situation has been already considered even though this model has only one stable 
particle. After setting the Boltzmann equations, we construct a simple model with two 
hidden DM sectors corresponding to two new f/(l) gauge symmetries, only one of which 
interacts with the SM sector. This is achieved through kinetic mixing of the new and 
standard gauge bosons. The DM particles are assumed to be Dirac fermions. We then show 
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how both particles can contribute to the rehc density of DM while only one can scatter 
elastically on nucleons. Although the direct detection rate tends to be rather high in this 
model, we show examples, where all constraints can be satisfied, including a case with a DM 
candidate around 10 GeV. 

This paper is organized as follows. The basic set-up is presented in section 2. An explicit 
model is constructed in section 3 and the implications for the relic density of DM as well as 
for the direct detection rates on nucleons are studied. Section 4 contains our conclusions. 



II. BASIC SET-UP 



We consider the case where xi and X2 are stable dark matter candidate particles. This 
can be achieved, for example, with a ® ^2 symmetry. We assume that m2 > nii with 
rrii = m^- and that X2 can only annihilate into Xi and not into SM particles. 

A set of coupled Boltzmann equations describe the evolution of the number density rij of 
particle Xi- I^i the following, X stands for some SM particles and rrix < frii. 
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where H is the Hubble parameter and n'^^ is the equilibrium number density of particle i. In 
solving the Boltzmann equations (II]) and ([2]), it is useful to introduce the variable Yi = ni/s 
describing the actual number of particle i per comoving volume, where s is the entropy 
density of the Universe. Solving these coupled Boltzmann equations, one can find Yi as a 
function of x = mi/T} In the new variables, the Boltzmann equations are recast as 
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^ Note that as opposed to the famihar case of coannihilation [15|, one cannot assume that n'^ jif^ ~ n\jni 
which rehes on the fact that scattering with background SM particles, e.g. X2X — >■ xi^'j occurs at a 
high rate. In the assisted freeze-out scenario, such reactions do not take place. Thus, we do not use this 
approximation. 
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where 



Hi 



{av)ij^ki{x) . (5) 

x=l 



The equihbrium number of particle i per comoving volume Y^^ = nf^/s has the following 
forms: 

Yr=^^yK,[x], Y,^^ = ^^^{rxrK,[rx], (6) 

g*s g*s '^T^ 

where r = m2/mi, Qi is the number of internal degrees of freedom of particle i and K2[x] is 
the modified Bessel function. 

The heavier particle X2 will not be in thermal equilibrium during freeze-out unless 
(at') 22^11 ^ {(yv)ii^xx- Thus, solving the Boltzmann equation in the case where many 
different processes contribute to annihilation would be complicated. In this analysis, we 
assume that s— wave annihilation processes dominate. This allows to simplify the compu- 
tation of the thermally averaged cross section for X2 when solving the Boltzmann equation 
while illustrating the dependence of the relic density on each parameter. 

If we limit our analysis to s-wave annihilation for simplicity, we can simply express the 
relevant s-wave matrix elements as 

a = yW22->ii = A^ii^22 , l3 = Mii^xx- (7) 

It is then straightforward to solve the two coupled Boltzmann equations ([3]) and ^ to obtain 
the abundance of Xi ^i-nd X2- To illustrate the assisted freeze-out mechanism, we choose the 
DM masses mi = 100 GeV and m2 = 150 GeV, taking different values for the annihilation 
amplitudes a and /3. The evolution of the abundances Yi, Y2 are displayed in fig. [T]for a and 
(5: (q!,/3) = (0.1,1), (0.01,1) and (1,0.1). The blue and red solid lines are respectively Yi 
and 12, and dashed lines show the corresponding equilibrium comoving number densities. 
When X2 interacts weakly with xi (see the top panels of fig. [I]), the abundance of xi is 
mainly determined by its interaction with SM particles and due to the weak interactions 
the final DM abundance is dominated with X2- A comparison of the cases a = 0.01 and 
a = 0.1 shows that the abundances of both particles are reduced by the interactions between 
Xi and X2- When the interactions of xi with X2 are stronger than with SM particles, i.e. 
a ^ (3 (the bottom panel of fig. [I]), the freeze-out of X2 is delayed and its abundance is, 
consequently, much reduced. Because of the interactions with X2, the abundance of xi 
increases in comparison with the previous case. The abundance of Xi now largely dominates 
over that of X2- 




FIG. 1: The evolution of the abundances of xi (blue) and X2 (red) per comoving volume as a 
function of x = nii/T. Dashed lines are the equilibrium comoving number densities. Dark matter 
masses are fixed as nii = 100 GeV and m2 = 150 GeV. The corresponding matrix elements are 
shown in each figure. 

III. DARK MATTER CANDIDATES WITH TWO EXTRA C/(1)'S 

As an explicit example of how the assisted freeze-out mechanism can occur in a specific 
DM model, we consider a toy model with a hidden sector containing two extra Abelian gauge 
symmetries, f/(l)' and U{1)", and two Dirac fermions ipi and ip2- The particle ipi is charged 
under both U{1)' and U{1)" gauge symmetries, and ip2 is only charged under the U{1)" 
symmetry. If ipi is the lightest particle charged under the U{iy symmetry and ip2 is the 
lightest one only charged under the U{1)" symmetry, then ipi and can be naturally stable 
particles like the electron in the SM. We assume that the hidden sector couples to the SM 
sector only through a kinetic mixing between U{1)' and f/(l)y. Then, the full Lagrangian 
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including this kinetic mixing is 

= ^SM - isine5,.X'-^ - Ix^J"^-^ - \x%X"'^'^ + \m%X'' + \m%,X"' 

- gx'X'^^il'^iJi - ^x"^^Vi7^V^i - 9x"X'l^^2'l''iJ2 + "^i^i^^i + "^2V^2^2 , (8) 

where the hidden f/(l)'s are assumed to be spontaneously broken leading to the gauge 
boson masses m-^, and m^„. In this toy model, we will neglect the possible effect of a 
mixing between the SM and hidden sector Higgs fields for the simplicity of analysis. In the 
SM sector, the mass of the Z gauge boson is and the gauge couplings are denoted by 
g = e/sy^ and g' = e/c^. 

The kinetic and mass mixing terms are diagonalized away by the following transformation: 

B = c^A - {t.Si: + s^c^)Z + {s^s^ - t,Ci.)Z' , 
= s^A + c^Ci:Z - c^s^Z' , 

X' = £iz + ^z' 

X" = Z", (9) 

where the rotation angle ^ is determined by 



'm%s,}r sin 2e 

tan 2e = 2 TT^ (1^) 



and the weak mixing angle s^r is very close to the physical value sw due to the stringent p 
parameter constraint. Then, the Z and Z' gauge bosons obtain the redefined masses: 

= + ^wMQ > (11) 



m 



2 



(12) 



on the other hand, the Z" gauge boson mass is the same as the X" gauge boson: m|„ = ''^i^,, • 
Let us list all the interaction vertices of the W, Z, Z' and Z" gauge bosons relevant for 
our analysis. The Z" boson has the following non-modified couplings: 

£ = -gx"Z';^[^p^J^^i,^ + i^^^H . (13) 

However, the other gauge bosons have modified couplings. In order to describe the interac- 
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tion vertices of W, Z and Z', let us define tfie various couplings, g^s, as follows: 
C = W+gJ[v^^PLe + u^PLd] + c.c. 

+ z^ [g% frPif + gfn h'PRf + gl ^in^i] + g^zw^W-]] 
+ [gf^ f^PLf + gf'n h'PRf + i^iYi^i] + gl^ [[z'w^w~]] 

+ [g^,, Z,Z^ + + g'k^zZ'.Z'^] ■ (14) 

One can find these redefined couplings expressed by the physical observables (unhatted 
parameters) in the appendix of 16| which describes a hidden sector with only a U{1)' sym- 
metry. The only difference between the gauge boson sectors of two models is the additional 
gauge boson associated with U{1)" that we assume to be decoupled from the SM sector. 

A. Relic density 

The relic abundances of two stable particles il)i^2 are determined by solving two coupled 
Boltzmann equations ([3]) and (jl]). The thermal relic density of DM is the sum of the relic 
densities of the two candidates ipi^2^ ^dm^^ = f^Vi^^ + ^^2^^- The free parameters of 
the model include the mass parameters mi,m2,mz' and mz" as well as the hidden gauge 
couplings gx' and gx" and the kinetic mixing parameter sine. In order to examine the 
dependence of the relic abundance of ipi^2 on each parameter, in our numerical anal ysis we 



search for the parameter space which satisfies the observed DM relic density limit 17|, in 
the mi — mz" plane fixing the other parameters.^ 

The annihilation of tpi pairs into SM particles proceeds through s— channel exchange of 
Z and Z': possible final states are //, W^W~ and Zh. The dominant annihilation mode 
depends on the DM mass mi, into only fermions when mi < M^y into additionally gauge 
bosons otherwise. In addition, all amplitudes are proportional to the kinetic mixing sine 
and the coupling gx'- The s— channel exchange can be strongly enhanced by a resonance 
effect. We therefore expect the relic density to drop rapidly when mi ^ m^/2 or m^//2. The 



^ In all the numerical analysis, we assume that mz" > 2toi to avoid a stable Z" . Moreover, for the Z' boson, 
electroweak precision test bounds and LHC detection prospects were already studied in Ref. HQ. We 



fix the kinetic mixing to its experimental upper bound given in 16l |18|. Smaller values oi gx' sine are, of 
course, even less constrained. 



7 



annihilation of X2 proceeds uniquely through the exchange of Z" and is therefore determined 
solely by the coupling gx" and the masses mz",m2. 

The relic densities l^^^/i^ and VL^^h"^ are displayed in fig. |21 in the mi —mz" plane for the 
reference case of gx' = 0.5, t^x" = 0.9, m<2 = 150 GeV and niz' = 150 GeV. The enhanced 
annihilation near the Z {Z') resonances explains the drop in Vt^^K^ for mi 45 (75) GeV. 
Similarly, the enhanced annihilation near the Z" resonance means that Vt^^h'^ is small for 
mzii ~ 300 GeV. Note that there is a lower bound, mz" > 2mi because we impose the 
condition that the Z" is not stable. The total DM relic density fi^j/i^ + fi^j/i^ is displayed in 
the bottom panel of fig. [2J The region between two thick dashed lines represents the points 
consistent with the recent WMAP relic density result [l3], allowing for 3a experimental 
uncertainties. In two bands at low values of mz" where mi ^ 35 GeV and mi ^110 GeV, 
the DM relic density is dominated by ipi while VL^^K^ dominates for mz" ~ 1 TeV. In the 
transition between these regions, both particles give a significant contribution to the total 
DM relic density. 

The annihilations of ipi and ?/'2 are controlled by the hidden gauge couplings gx' and gx", 
respectively. To study their effect, we reduce each coupling separately. In fig. [3l the total 
DM thermal relic density fi^^/i^ + ^^^2^^ is shown for the two representative cases gx' = 0.3 
or gx" = 0.3 in the mi — m^// parameter space. For the other parameters, we take the 
same values as in the reference case of fig. |2l Comparing the bottom panel of fig. [2] with 
{gx'iQx") = (0.5,0.9) with the left panel of fig. [3] with {gx',gx") = (0.3,0.9), one sees that 
the allowed bands at low values of rriz" move closer to either mz/2 or mz'/2. This is because 
as the hidden gauge coupling gx' decreases, the annihilation of tpi is weaker, and stronger 
resonance effect is therefore required to obtain appropriate annihilation strength. This also 
means that in the mass region between the two resonances, mi ^ 60 GeV, the relic density 
of ipi increases. Therefore, the total DM density can be in agreement with the measured 
value for lighter Z" masses. When the hidden gauge coupling is decreased to gx" = 0.3, 
the requirement of a stronger resonance effect means that the region in agreement with the 
observed value of the relic density moves to lower values of rriz", see the right panel of fig. |3l 

The interactions of the DM particle tpi also depend on the Z' mass m^/. In fig. HI we 
display the total DM relic abundance ^^^h"^ + ^^2^"^ "^i " ''^z" parameter space to 

illustrate the dependence on the Z' mass m^/. In the figure, the left panel corresponds 
to nizi = 200 GeV and the right panel to m^/ = 25 GeV. In this analysis, we used the 



8 




FIG. 2: Contour plots for the relic abundances of the dark matter particles ipi (top-left), "02 (top- 
right) and the total (bottom) in the mi — mz" plane. We fix the other parameters as follows: 
9x' = 0.5, gx" = 0.9, m2 = 150 GeV and mz' = 150 GeV. In the bottom panel, the region 
between two thick dashed lines is allowed by the WMAP result on the DM relic density. 

reference values for the other parameters. As can be seen clearly from fig. |U there are two 
well separated resonance regions, one around mz/2, the other near mz'/2. In addition, the 
right panel shows that the lighter particle ipi can be a DM candidate around 10 GeV when 
Z' is light. Thus, this model could have a light DM candidate as hinted by some direct 
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FIG. 3: Contour plots for the total relic abundance of the dark matter particles ipi and "02 in 
the mi — mz" plane. The left and right panels correspond to {gx',gx") = (0.3,0.9) and (0.5,0.3) 
respectively. The other parameters are fixed at the reference parameter values. In each plane, the 
region between two thick dashed lines is preferred by the WMAP DM relic density result. 

detection results, see more details in the next section. 

Finally, we investigate the dependence of the DM relic density on the ip2 niass In 
fig. [5l the total DM relic density of ipi and ip2 is shown for the representative case m2 = 200 
GeV in the rrii —mz» plane which is to be compared with fig. [2] for which m2 = 150 GeV. All 
other parameters are kept to the reference values, gx' = 0.5, gx" = 0.9 and m^/ = 150 GeV. 
As expected, for larger 1^2 niass, the preferred parameter space moves to larger mz" region, 
now mz" ~ 1.2 TeV due to the change of the mz" resonance region while the preferred value 
for mi at low values of m^" is not much shifted. One can see the tendency from the bottom 
panel of fig. |2] and fig. [5l 

B. Direct detection 

We have shown that this simple model realizing the assisted freeze-out mechanism can 
satisfy the observed DM relic abundance. We now discuss the influence of this mechanism 
on the prospects of observing DM in direct detection experiments. The most distinguishing 
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FIG. 4: Contour plots for the total relic abundance of the dark matter particles -01 and 11:2 in the 
nil — rnz" plane for the cases mz' = 200 GeV and rriz' = 25 GeV. The parameters gx' , Qx" and 
m2 are fixed as the reference values. The regions between two thick dashed lines are allowed by 
the recent DM relic density observation. 

feature of the assisted freeze-out scenario compared with the other multi-DM scenarios is 
that only lighter DM particle ipi can directly interact with the SM sector. Consequently, 
only can be detected in DM direct detection experiments. The lighter DM particle 
can elastically scatter off a target nucleus through t— channel Z and Z' gauge boson 
exchange. One can easily calculate the spin-independent (SI) ^/'i-nucleon cross section using 
the following effective operator: 

A// = ()/*i7^'i*i/7''/, (15) 

where = sil^il^M + ^IMt±»Sl . (16) 

^ 2m| 2m|, ^ ' 

The current experimental bounds on are extracted from DM direct detection ex- 
periment results assuming that the couplings to protons (/p) and neutrons (/„) are equal. 
However, the couplings are different in this model. In order to compare directly with the 
limits on given by experiments, we thus use the normalized cross section on a point-like 
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FIG. 5: Contour plots for the total relic abundance of the dark matter particles il^i and ^2 in the 
mi — mz" plane for the case m2 = 200 GeV. The reference values are used for the other parameters. 
The regions between two thick dashed lines are allowed by the recent DM relic density observation. 



nucleus |l9| : 



SI /i^, [Zf, + {A-Z)l^]' 

Moreover, the experimental limits are also extracted assuming that the local DM density is 
due to only one DM candidate. However, in our model only the lighter particle ipi can be 
observed by DM direct detection experiments. Therefore, we rescale the SI scattering cross 
section by Q^t-^h"^ /QBuh"^ assuming that the contribution of each DM particles to the local 
density is the same as their contribution to the relic density. 

In fig. ini we present the normalized SI scattering cross sections of tpi as a function of rrii 
for (gx', mz>) = (0.5, 150 GeV), (0.5, 200 GeV) and (0.3, 200 GeV) fixing gx" = 0.9 and 
m2 = 150 GeV. The value of mz" is chosen to follow the contour VL^-^h? + VL^^h? = 0.13. 
The first feature is that the rescaled cross section decreases around the resonance points 
(mi ^ mz/2, mzi/2), this is simply because VL^p^h? decreases sharply near the resonances. 
Actually, the elastic cross section itself has a much milder and smoother dependence on mi. 

Comparing the cases {gx'^f^z') = (0.5, 150 GeV) and (0.5, 200 GeV), one can see that the 
cross sections become smaller for larger mz' because of the suppressed contribution of the Z' 
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FIG. 6: SI scattering cross section of V'l normalized to Xenon target nucleus and rescaled by 
O^^/i^/r^DM^^ (see text). The red solid, green dotted and blue dot-dashed lines correspond to 
{gx', mz') = (0.5, 150 GeV), (0.5, 200 GeV) and (0.3, 200 GeV), respectively, gx" and m2 are 



20(], is shown 



fixed to 0.9 and 150 GeV. The experimental limit, which is taken from XENONIOO 
by the black dashed line. 

exchange diagram. Thus, one can recover a region where the elastic cross section drops below 
the limit of XENONIOO 20j. Both the Z and Z' exchange diagrams are also suppressed for 
smaller values of gx', we therefore obtain smaller cross sections in this case as can be seen 
by comparing the curves for {gx'.'^nz') = (0.5, 200 GeV) and (0.3, 200 GeV). Note that in 
fig. [6] there is a discontinuity in the direct detection curves when rriz' = 200 GeV, which is 
because the Qh"^ = 0.13 contour is also disconnected, see fig. HI In our analysis, we fix sine 
to the experimental upper limit as mentioned in the previous subsection. If we use smaller 
sine, the effect is very similar to the case of smaller gx' since the interactions between ipi 
and SM particles are approximately proportional to gx' sin e. In summary, we therefore find 
that it is easier to satisfy the direct detection constraint for smaller values of gx' sin e and 
larger values of mz' than those of the reference case. 

In this analysis, we apply the normalization of the cross section to Xenon, Eq. (fT7|) . since 
XENONIOO provides the most stringent limit in most mass range. Actually, the cross section 
for Germanium is very similar to the Xenon case due to similar Z/A ratios in both nuclei. 
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FIG. 7: SI scattering cross section of tpi for the Xenon target with gx' = 0.5, gx" = 0.9, m2 



150 GeV and m^' = 25 GeV. The experimental Hmit from XENONIOO 
dashed hne. 



20| is shown by the black 



In order to check the possibility of a light DM hinted by the data from DAMA 10| . 



CoGeNT 



ll| and CRESST [l^, we show the SI scattering cross section corresponding to 



a light Z'. We choose niz' = 25 GeV using gx' = 0.5, gx" = 0.9 and m2 = 150 GeV in 
fig. [71 As can be seen from the figure, the scattering cross section is severely constrained 
by the XENONIOO limit since the scattering rate is enhanced due to the small mass of Z'. 
Nevertheless, we can find small allowed region around the Z' resonance point of mi ^ rriz'/^. 
One can easily expand this region, where the SI cross section drops below the XENONIOO 
bound, using smaller gx' sin e. In addition, lighter DM particle ipi can satisfy the scattering 
cross sections required by the DAMA, CoGeNT or CRESST results. 



IV. CONCLUSION 



We have illustrated with a simple toy model containing two stable dark matter particles, 
how the assisted freeze-out mechanism worked and could reproduce the measured value for 
the relic density of dark matter. The main feature of this type of model is that only one 
of the DM particles is involved in direct detection searches while both contribute to the 
relic density. In particular, when the DM particle that interacts with SM particles is the 
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subdominant DM component, it is possible to reconcile models with large elastic scattering 
rates on nuclei with the exclusion bounds of XENONIOO. Moreover, the lighter particle can 
be a light DM candidate around 10 GeV as indicated by the DAMA, CoGeNT and CRESST 
results. 
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